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ABSTRACT 



A Carousel Remotely Piloted Vehicle (RPV) Recovery 
System consists of a vertical mast which, supports a horizon 
tal member that is free to rotate. The approaching RPV is 
caught by vertical cables suspended from the horizontal 
member and the kinetic energy of the RPV is dissipated 
through the motion of the recovery system. This thesis 
presents a simplified dynamic analysis to describe the 
motion of an RPV after impact v/ith a recovery system simu- 
lating the carousel system. The eo^uations of motion for 
the RPV were obtained from Lagrange’s equations and modeled 
using the _ Continuous System Modeling Program on the IBM 
360-67 computer at the W.R. Church Computer Center, BPS. 

The results showed that damping was required on motion of 
both the horizontal member and the suspended cables in orde 
to prevent possible damage to the 150 pound RPV at an as- 
sumed 50 knot approach speed. A sensitivity analysis was 
performed by varying the system design parameters and the 
results are presented in the body of the thesis. 
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I. INTRODUCTION 



Extensive efforts on Remotely Piloted Vehicles (RPV) 
have been underway for many years, covering a v/ide range of 
applications and missions, Ref. 1. It is only recently, 
however, that both the need and the technological capability 
for mini-RPVs have coalesced to the point that the develop- 
ment of these systems has become imperative. A mini-RPV 
is generally defined as an RPV that weighs less than 200 
lbs mission loaded. 

An essential part of an RPV system is the recovery 
system. The Army has recently experienced several mini- 
RPV losses due either directly or indirectly -to the 'recovery 
method. The methods used to recover mini-RPVs vary from 
catching the RPV with an arresting cable similar in prin- 
ciple to that used on an aircraft carrier, to guiding the 
mini-RPV directly into a vertical net which is designed to 
flex and absorb the kinetic energy of the RPV. A completely 
satisfactory recovery system has not yet been established, 
hence studies into alternate recovery systems are of con- 
tinuing interest to the Army. 

The purpose of the analyses described herein wa.s to 
model the dynamics of one alternate recovery system candi- 
date, the "caroLisel" method. The carousel method was 
initially described as a candidate for a recovery procedure 
in a concept stiidy performed by Teledyne Ryan for the Army, 
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Ref. 2, but the initial concept formulation did not include 
significant dynamic considerations. Therefore, the analyses 
conducted in this thesis, which are primarily oriented to- 
wards system dynamic response behavior, should serve in the 
role of providing additional information for concept feasi- 
bility considerations. 

The carousel recovery system, shown in Figure 1, was 
conceptualized in this analysis as being composed of a 
vertical mast that supported a horizontal member which v;as 
allowed to rotate about the vertical mast. Hanging from 
the horizontal member are cables that will be caught by 
the wing tips of the approaching RPV. The RPV’s kinetic 
energy will be dissipated by the recovery system through 
various me'ans that will be evaluated in this analysis. 

The RPV must be brought to rest without allov/ing 
excessive RPV swing since excessive swinging motion could 
result in the RPV hitting the recovery system structure. 

This consideration was of paramoimt concern when evaluating 
the merits of recovery system parameter variations. 

For the purpose of this analysis the mini-RPV to be 
recovered was assixmed to be the AQUILA, which is a U.3. 

Army mini-RPV developed by Lockheed liissiles and Space 
Company, Ref. 1, and is sketched in Figure 2. It has the 
follov/ing nominal characteristics: v/ing span, 12 ft; 

weight 150 lb; approach speed 50 knots; maximum loading, 

6g vertical and axial, 3g lateral. 
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Figure 1 Carousel Recovery System 
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Figure 2 Aquila Mini-RPV 



II. ANALYSIS 



A. ASSUIvrPIIONS 

Simplifying assiimptions, as defined below, were made 
when establishing the equations of motion using the pro- 
cedures suggested by Lagrange's equations, Ref. 3> in order 
to establish the basic principles of the carousel recovery 
system dynamics. Inclusion of system details such as re- 
covery boom elasticity and RPV mass moments of inertia were 
considered as items which might tend to mask the basic 
question of feasibility at the present stage of analysis. 
Under these assumptions, the results supported the viability 
of the carousel recovery method and a basis was established 
upon which further analyses may be performed when the sim- 
plifying assumptions are removed. 

The recovery system was initially assumed to be con- 
servative, that is, friction losses were ignored and as 
soon as the RPV contacted the recovery system the engine 
was secured. This assumption was subsequently changed by 
the addition of damping to the recovery system configuration. 
The horizontal member wa,s assumed to be rigid. A careful 
choice of material and structural characteristics for the 
horizontal member can minimize its dynamics and allow it to 
be approximated by a rigid member. The cables hanging from 
the horizontal member v;ere replaxed by massless thin re- 
covery rods. The RPV v/as treated as a point mass and v;a,s 
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captured by one recovery rod. Treating the RPV as a point 
mass simplified the coordinate system used in the analysis 
by allowing the mass moment of inertia of the RPV about its 
principal axes to be neglected. The RPV v/as also assumed to 
approach perpendicular to the plane formed by the horizon- 
tal member and the vertical recovery rod that cauight it. 

The approach velocity of the RPV wa,s assumed to transfer 
undiminished to the recovery rod giving the rod an initial 
angular velocity only in the 0 direction. Reference should 
be made to Figure 3 v/here the orientation of the generalised 
coordinates 0, 0 and Ip' is portrayed. These coordinates 
are a form of Euler angles, independent of each other, that 
together completely specify the position of the RPV relative 
to the fixed reference frame x-y-z. First 0 defines the 
rotation of the moving coordinate system x'-y'-z' aboiit the 
fixed axis z. It describes the rotation of the horizontal 
member in the x-y plane. Next 0 defines the motion of the 
RPV about x' v/ithin the moving coordinate system. This is 
portrayed in Figure 3 as the motion of the RPV projected on 
the y'-z' plane. Finally, Ip" depicts the remaining motion 
of the RPV after it rotates through 0 and 0 . 

The system identified in Figure 3 is the reference or 
base configuration that was assumed for the analysis based 
on judgement of the RPV's weight and size. In this refer- 
ence conf iguration the horizontal member, R^, was 28 ft 

long and made of 6-inch 8,luminum I-beam having a mass 

2 

moment of inertia of 1159 ft-lb-sec about pivot "0". 
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Figure 3 Simplified Recovery System 
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The recovery rod R^, was 6 ft long. The RPV of mass 4.66 



slugs (150 lbs) was suspended at the lower end of R^. V/ith 

an assumed approach speed of 50 knots (84.4 ft/sec), the 

♦ 

initial angular velocity, 0^, of the massless rod R 2 was 
14 rad/sec according to the relation of equation (l) : 



A linear, torsional spring, of spring constant K, lo- 
cated at "O'*, the pivot of the horizontal member, v/as 
initially used to resist the rotation of the horizontal 
member, v/ithout any damping in this reference configura- 
tion it v/as further assumed that the spring's stored energy 
could be eliminated from the system through a ratchet re- 
lease mechanism, after the horizontal member had swung 
through its initial maximum displacement. The ratchet 
mechanism was not modeled, as it was felt, and later con- 
firmed by the results, tha,t the critical motion of the RPV 
would occur shortly after impact with the recovery system. 
Simulation of the ratchet mechanism would have no effect 
on these results. 

B. OVERVIEW OP AIhA.LYSIS 

The general approach of the analysis was to model the 
RPV motion after engagement with the recovery system using 
the generalized coordinate system Q-Cp-lp' indicated in 
Pig. 3 . If fbe assumptions of a rigid horizontal member 
and a point mass RPV representation were removed, a minimum 




( 1 ) 
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of five more generalised coordinates would have been re- 



quired. The equations of motion for the three degree of 
freedom system were obtained using Lagrange's equations 
for a conservative system. From these equations the 
a.cceleration, velocity and position of the RPV as a function 
of time v/ere obtained in component form along each coordi- 
nate direction through the use of the Continuous System 
Modeling Program (CSivlP) . The input parameters were then 
varied and their effect on resulting dynamic motion studied. 

Changes to the basic recovery system which included the 
addition of damping terms were later entertained as a con- 
sequence of initial results. 



C. DETAILED ANALYSIS 

1. Lagrange ' s Equation 

The equations of motion were established using 
Lagra,nge's equation, Ref. 3. Three ec^uations of motion 
were obtained for the system, since there were by assi-unp- 
tion only three generalized coordinates: the coordinates 

0 , 0 and t/y ^ as defined in Fig. 3. Lagrange ' s equation 
in the conservative form, without any generalized forces, is 



-J: ( Al 1 _ Al 

dt y^qj_ j ^q^_ 



i = 1, 2, 3, 



( 2 ) 



where 

q^ = generalized coordinate Q ,0^=0 

Q. = time derivative of q. coordinate 
-1 
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T = T(qj_,4j_) = Kinetic energy (ft-lb) 
V = V (cij^) = Potential energy (ft-lb) 



The kinetic energy was defined rela,tive to the inertia,! 
reference frame by the following expression: 



T 







( 3 ) 



where 

I 

o 

I.I 

e 



V 



SPV 



= Moment of Inertia of the horizontal member 

2 

a,bout pivot point "0” (ft-lb-sec ) 

= Mass of HPV (slugs) 

= Angular velocity in 0 coordinate direction 
(rad/sec) 

= Velocity expression for RPV derived in the 
next section (ft/sec) 



The potentia.1 energy included the gravitational term acting 
upon the RPV mass and the stored energy in the torsional 
constra,int spring. Equation 4 expresses the potential 
energy relative to the inertial frame of reference as 
follows: 

V = MgR2(l-cos0cosT/r^+iK0^ (4) 

v.'here 

R2 = length of recovery rod corresponding to RPV 

distance from the horizontal member follow'ing 
recovery, (ft) 

K = spring constant of spring at pivot of horizontal 
member (ft-lb) 
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Since its exactness is critical to tlie evaluation of 
Lagrange's equations of motion, the RPV velocity is treated 
separately in the ensuing section using the generalized 
coordinates . 

2 , Derivation of RPV Velocity 

To derive an expression for the RPV's velocity 
after impact with the recovery system, the position of the 
RPV within the moving x*-y'-z' coordinate system as a func- 
tion of the generalized coordinates 0, l/f was transformed 
to the fixed inertial reference system x-y-z. The time 
derivative of the RPV's position, expressed as a function 
of the three generalized coordinates 0,0 and 0^ and ref- 
erenced to the fixed inertial system, represented the ah- 
•solute' velocity of the RPV. 

The RPV's position within the moving frame x'-y'-z' 
was expressed as follows: 



x’ = -R 2 siil 


( 5a) 


y' = R 2 sin0cos0‘ 


(5b) 


z' = -R 2 COS0COS0" 


(5c) 



Transformation to the fixed inertial coordinate 
system may be stated as: 



r b 




— 




1^- 


> = 


COS0 


- sin0 


y 




sin0 


cos 0 






— 


— 


z 


r 


s' = 


-R^CO£0 



r 



-RpSin ’\j/~ 






R^COS0 



+ < 



R2Sin0cos0~ R-,sin0 

J L “ J 



(6a) 

(6b) 
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In an alternate form, the three components expressing the 
absolute position of the RPV mass following recovery en- 
gagement may be expressed as: 



X 


= RpCos 0 -R 2 sin"^cos 0 -R 2 Sin 0 cosl/Osin 0 


(7a) 


y 


= R^cos0 -R^s±n1p~sinQ +R^s±n0 cos'lp'cos 0 


(7b) 


z 


= -R^cosCpcos ilf 


(7c) 


The absolute velocity of the PlPV was obtained by taking the 


first 


time derivative of equation 7 , thereby yielding 


# 

• 


X 


= j^-R^sin 0 +R 2 ( sinl/^sin 0 -siii 0 cos 0 ’cos 0 )J 0 






-(R^cos^ cos Ip'slnQ )(p +R 2 ( sin 0 sinl/fsin 0 - 






cosT/Ocos 0 ) ijf 


( 8 a) 


y 


= ■ ^R^cos 0 -R^C sin 0 *cos 0 +sin 0 cosl|/rsin 0 ^ 0 






+ (R^cos0cos1//“cos 0 )(p -R^ ^cos 0“sin 0 +sin0 


slnlj/' 




cos 8 ^ijr 


( 8 b)' 


Z 


= Rp [(sin 0 cos 0 ') 0 +(cos 0 sin 0 ') 


( 8 c) 




The expression desired for use in Lagrange's 
9 





2 

eqLiation, obtained by squaring the components 

separately then adding them; i.e., 



„ 2 .2 .2 .2 

^RPV = - +y +z (9) 

• 2 -2 -2 

The expressions for x , y and z are equations 10a, b 
and c . 
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= |^-R^sin 0 +R 2 ( sin'l/z'sin 0 -sin 0 cosl^cos 0 )J ^ 0 ^ 

o ^ — 0*0 

h-R2^cos ^cos^Y^sin*" 0 0 

+R 2 ^^( sin 0 sin 0 'sin 0 -cos T/^cos 0 ) 1 /T 

- 2 R 2 COS 0 cos 0 'sin 0 ^-R^sin 0 ^-R^C sin'^sin 0 -sin 0 

coslIfcosQ )] 0 0 

+ 2R^{s±n0s±nlj/'s±nQ -cos'\l/'cos Q ) |^-R^ sin 0 +R^( sinyr 
sin 0 -sin 0 cosTj/tcos 0 ) 0 ijf 

- 2 R 2 ^cos 0 cos' 0 'sin 0 ( sin 0 sinl/^sin 0 -oos'lp'cosQ )<p')j/' 

(10a) 



.2 

TT 

c / 



|r^cos 0 -R^C sinl//*cos0 +sin0cos0'sin 0 )J-^ 

2 2 2 ' 2 
cos 0COS l/rcoG 0 0 

p r n 0 * 0 

+R 2 cosl^sin 0 +sin0 sinT^/rcos 0 I 0^ 

+ 2 R 2 COS 0cos0'coG 0 jR-| COG0 -R^C sinl/Tcos 0 +sin 0 
cos'l/^sin0 )]0 0 

- 2 R 2 jcoG0'sin0 +sin0sin0'cos 0j j^R^cos 0 -R 2 ( sinl/r 

COS0 +sin0cos0'sin 0 J 00’ 

2 

-2R2 c os0cos0'cos0 jjcos0'sin0 +sin0sin0' 

(10b) 
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2r2 2 2 2 2*2 

pin 0 cos llr(p +COS 0 sin 0- +2sin0cos0 



sin T/^cos"^ 0 



(10c) 



• 2 -2 *2 

When X , y , z v;ere combined, several cancellations 

2 

occurred and the final expression for simplified to; 

XlX V 



V 



RPV 



j^R^^+R^^ sin^0"+R2^sin^0 cos^0'-2R^R2Sin0'j 0^ 
+R2^cos20'0^+R2^1^+2R2Cos0 cos •^'(R^-R^ sin 0*) 0 0 



2 . 



+ (-2R]_^2^^^0 sin0)00~ 



( 11 ) 



3 . Derivation of Equations of Motion 
a. 0 Equation of Motion 

2 

Now that an expression for Vj^py has been 
obtained, Lagrange's equation can be evaluated according to 
section II.G.l. The equation was evaluated for each gener- 
alized coordinate 0 , 0 and ijr , The 0 equation of motion 
was obtained by evaluating the Lagrange's equation using 
the 0 generalized coordinate. 



d rhT 






30 



£2 + II = 0 

30 W 



( 12 ) 



A^^ain: 



Kinetic energy, T=lrlQ 0 ^+-|-IvIVppY^ 

Potential energy, V=i5gR2 [l-cos0 cos +-J-K 0 



(13) 

(14) 



1=0 

30 
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(15) 



( 16 ) 



+M [(R^^+R2^sin‘^-^+R2^sin^0cos^1^-2R^R2sinYr)0 



2 . 2 



2 . 2 



+R2Cos0cosT/^(R^-R 2 sin'|//')0 +(-R^R2Sin0 sinl^+ 

2 • <p)iir 



r,2 sin 



2 * 2 2 * 
+ (2R2 sinYrcosl^7^+2R2 sin0cos0cos 'ijf 0 

-2R2 sin 0 coslZ/^sinl^ T^-2R^R2 Cos-j^ 1^) 0 



-R2sin0cosl/r(R^-R2sinl^) 0 -R^slnllfcos 0{Rj^-R^ 

sinl/^) 0 ij/ 

2 * • • . . 

-R2 cos0'cos0cos0'0'0-RiR2Cos 0 sin"^0 0'-R]_R2 

* 2 

sin0cos0'0' 



+R2"^cos0 01^ 



( 17 ) 



Substituting eouiations 15 , I6 and 17 into equation 12 and 

2 

dividing through by (I.IR2 ) yields the following 0equation 
of motion; 
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?.© 



Z R 



R. 

+s±n(p { 1-^3 ±n1p') l/r 
Rn 

+ COS0 (^cosi//'--^-sin2T/r) 0 

2 R. 

+ sin2T//~cos 0 -2^cosl/r) 0 



* * 

+ sin20 cos ■}p" 0 q 



+COS0 (-2^sim//'-cos20'+l) 0 0- 

R 

1 • 2 
-sin0 (tj— cos0"-vsin20‘) 0 



R 






; 2 



- •^siii0cosi/r0'' 



+ fW 2 = 0 



( 18 ) 



b, 0 Equation of Motion 

Similarly the 0 equation of motion was ob- 
ta,ined from the following Lagrange's equation: 



d / 3T 



- ll + ^ = 0 
20 20 



(19) 



where 
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d/5T 



dt 



h(p 



LI [^2^003^1^0 +R2COS0 cos'l^(R^-R2sinl/r) 0 

p • • • • 

-2R2^co B'ip' sinT^'00 -R2sin0cos l/^(R^-R2sin0')0 0 

• • _ * « 

-R2sinl/rcos0 (R^-R2sinT/r)00' -R^^oosCp cos2-^i^0j 

( 20 ) 



= M [R2^sin0cos0 cos'^-0-0'^-R2Sin0 cosl/r(R^-R2sini^)00 

2 * * 

+ (-R2R2COS0 sin1^+R2 cos0)01|//'j (21) 



.2., '2 



90 



LlgR 2 s in 0 c o s 



(22) 



Substituting equations 20, .21 and 22 into equation 19 and 

2 

dividing through by (IVIR 2 ) yields the following 0 eq_uation 
of motion: 



R- 



^cos0 coBlff-f-cos0 sin2'0' 



e 



+ COS 



't<p 



-sin2 0" l/r 0 

• • 

-COS0 ( CO s2 1 / 0 + 1 ) 0 i/r 

--o-sin20 cos^l/o 0^ 



+ ^sin0cosl/o= 0 



(23) 
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c. lf/‘ Equation of Motion 

The lp~ equation of motion was similarly obtained 
using the following Lagrange's equation: 



where 



(2l \ - + IX ^ 0 

\^f) dt 




(24) 



^ f 



= Ijlf+i- ^cos 0 sini^+cos 0 ) 00 - ^cos 0 'sin 00'0 

Rl 

+ (- ^sin 0 sin 1 ^+sin 0 ) 0 J 



(25) 



~ = R 2 ^M |^( sin 0 *cos 0 '-sin ^0 cosT^sinT^- •^cos 0 “) 0 ^ 

Or 2 

* 2 T 2 

-cosT/rsinT/r 0 +(- ^cos 0 sin 0 '(R^-R 2 sin'l/r) -cos ij/- 

COS0 ) 0 0 



R 



_ ^sin0cos0- 0 T^] 



( 26 ) 



^ = ^cos0s±n~\lr (27) 

Substituting equations 25, 26 and 27 into ea_uation 24 and 

2 

dividing through by the (lvIR 2 ) .yields the following 
equation of motion: 



23 



R 



(- ^sin0sinl^+sin0 )0 + 



+ COS0 (l+cos2'^) 0 0 



+ (^sin2 ■]//"( sin^0 -1)+ —cos'ij/') Q‘ 

2 

*0 rr 

+-^sin2l/'' 0^^+ ^cos0sin0'= 0 



R. 



( 28 ) 



4 . Solving the Bquations of Motion 

Th.e equations of motion were solved for time his- 
tory information of the RPV's position, velocity and accel- 
eration using the IBLI Continuous System Modeling Program 
(GSIIP) which was available as resident softv;are on the IBM 
360-67 computer located in the VV.R. Church Computer Center, 
Naval Postgraduate School. The use of CSLIP as a high-level 
programming tool is quite similar to the familiar engineer- 
ing language of PORTRArl, and. programming details may be found 
in the test by Speckhart and Green, Ref. 5. Details of the 
actual program used are presented in the Appendix. To use 
this software package, the three equations of motion were 

first solved for the highest order derivative 0 , 0 and if' 
in terms of the three generalized coordinates and their first 
time derivatives. Then, in a manner similar to an analog 
computer solution of differential equations, the CSJIP program 

proceeded to integrate these equations numerically for values 
for 0 , 0 , , 0,0 and lp~, which were fed back into the 

eo^uations for 0 , 0 and ijr . This itera,tive method of , 
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solution was repeated at ea.ch time increment. The resulting 
values for 0 , 0 and were listed as a function of time 
and from these the motion the RPV v;as studied to assess the 
feasibility of the recovery system when a particular set 
of structural parameters and initial RPV conditions was 
specified. 

Since the governing equations of motion (eqns. 18, 

23 and 28) were inertially cross coupled, it was necessary 



to rearrange the equations into an inertially uncoupled form 
that would be acceptable for solution by the CStlP. In matrix 
form, equations l8, 23 and 28 may be expressed as; 





















ail 


ai2 


ai3 




0 " 




H 




^21 


0 


^23 


< 


t 


> = < 


^2 ^ 




^“31 


^32 


0 




0 






whe ] 


re 
















I 


R 




, 






ail 


0 




p 

-slnljf) +CO 


2 2 
s ijfsln 0 





Rn 



^12 = 2-31 = siu0(l- ^sini^) 



R. 



‘'13 “ ^21 “ cos0(^cos0'-lrsin^l/d) 



(29) 



(30a) 

(30b) 

(30c) 



^23 = cos ijr 



a .^2 - ^ 



(30d) 

(30e) 



30 



(30f) 



R-, . p R 

= •^sin0cosl//'T//'^^+sin0 (^cosl//'-o-Gin2 l/r) 0 



Ri 



-(sin T/^cos 0 -2^cosT//')0i/t -sin20cos llf Q 0 
R-j • • F o 

-(l-cos2l/^-2^sin0-)cos00-0- 0 

b2 = ijcos^'0'sin2 0 0^+cos0 (1+COS2 0') 00 +sin20'1j^0- 



^sin 0 cosTp~ 



R. 



= (-^sin 2 0 'cos ‘^0 - ■^coslp') 0 ^-- 2 -sin 2 0*0 



cos 2 l /^)00 - ^cos 0 sin 0 ‘ 



(30g) 
-cos 0 ( 1 + 

(30h) 



Substituting equations 30a through 30h in matrix expression 
29» 0> 0 and ijf were solved using Cramers' rule, Ref. 4. 



( 31 a) 






T^= 1^1 0 = Ol 

^ ‘ A| (31b) ^ I A| (31c) 



wnere 



A 



^11 "'12 ^3 
3'21 ^ 

^31 ^32 



c^23 

0 



= ^12 ^23 ^31-^11 ^23’^^21 ^13 



( 32 a) 



A- 



°1 

bo 



"^12 13 

0 ag3 

1 0 



b^ a^2 ^23”^1 ^23'‘’^2 ^13 



( 32 b) 
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A. 



a 



a 



a^ 



a 



a 



11 




^13 






21 


°2 


a23 




^' 11^23 ®' 21 ®' 13 ^ 


31 


^3 


0 




( 32 c ) 


11 


^•12 


^1 






■21 


0 


^2 




a 2 i )+ b 3 


'31 


1 


^3 


(-a2iai2) 


( 32 d ) 



Squations 30a through. 32d v.^ere prograramed on the IBIvI 300-67 
computer using CSI.IP, the details of which are described in 
the Appendix, The program evaluated S , 0 , , Q 

0, 0 and' as a function of tine from the initial con- 
dition of RP'V impact. The evaluation of the computer 
generated time history results is described in the next 
section. 
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III. SYSTEM DESIGN ANALYSIS AND RESULTS 



A. UNDAI,rPED SYSTEM 

The reference or base configuration of the recovery 
system described in section II. A. and Pig. 3 v;as the first 
configuration considered in the analysis. Prom this con- 
figuration each parameter v;as varied individually for various 
RPV approach speeds to assess the system's sensitivity to 
their changes. 

Prom the outset of the investigation it v/as felt, as 
supported by early results, that the motion of the RPV 
along the 0 coordinate would be a primary system constraint. 
A guideline value for 0maximum ( 0max) was chosen to be 90 
degrees, v/hich corresponds to the recovery rod (R 2 ) reaching 
the horizontal. Por the recovery system to pa,ss its first 
requirement, an RPV approaching at 50 knots would have to 
be recovered v/ith the maximum amplitude of 0 not exceeding 
90 degrees. 

The results of modeling the reference configuration 
suggested that energy disipation would hsive to be included 
to bring 0max down to 90 degrees or less. V/ithout any 
damping the system exhibited unacceptably high magnitudes 
for 0 that corresponded to the RPV completely rotating 

e.bout the horizontal member when approach speeds exceeded 
25 knots. This can be seen by referring to Table I where 
0 is listed for the base configuration as a function of 
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VARIATION FOR UNDAMPED SYSTEM 
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time after impact for approa.ch speeds ranging from 7 Imots 



impact v/as greater than l80 degrees. The SPY has started 
to rotate over the Iiorizontc,! member. Va.riation of the 
system param.eters did not correct this deficiency''. 



In an effort to correct the excessive RPV sv;ing in the (p 

coordinate direction, a Rayleigh dissipative energy term, 

D, v;as introduced into Lagrange's equation: 



This e,dditional term precipita,ted through the equation 

of motion 6.erivation of section II. C. v;ith the addition of 

» 

only a Oq0 term to the 0 ec[uation of motion, equation 18. 
This cha,nge v;as easily incorpora,ted into the computer pro- 
gram as indicated in the Appendix. 

An initial estimate of the magnitude of the dam.ping 
constant v;as obta.ined using the simplified linearised 
pendulum equation belov;. The damping value v;as varied to 
detect the system's sensitivity to C^. 



( 0q= 2 rad/sec Ref. Equation 1) to 50 knots ( 0^=14 rad/sec). 
Note that for 0^=7 rad/sec, or approximately a 25-knot 
approach speed, the magnitude of 0 at 0.7 seconds after 





( 33 ) 



• 2 

where D=-lC^ 0 



0 +2 0 -r (SJ^ 0=0 



(34) 
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v/here 



0 

For values of 1^=1159 Ib-ft-sec^, K=50 ft-lb and <^= 0.15, 
v/iiich v/as selected to introduce just a sms,ll amount of 
da.mping, the initial estimate of "'^as 75 ft-lb-sec. 

She results of adding damping to the pivot of the 
horizontal member in the base conf iguration indicated that 

it v;s.s only partially successful in reducing 0man. This 
can be seen by referring to fable II where 0 is listed for 
this configuration as a function of time after impact for 
a.pproach speeds ranging from 7 knots, v;here 0q=2 rad/sec, 
to 50 knois, where 0^=14 rad/sec. Ilote that the magnitude 
of 0 was still approaching 180 degrees for approach speed 
initial conditions in excess of 35 knots ( 0^=10 rad/sec). 

Variation of s^/'stem parameters including did little 
to further reduce 0max, as shown in Figures 4 through 8 , 

where the base configuration of R-|l =28 ft, R 2=6 ft, LI= 4.66 
slugs, 1^=1159 ft-lbs-sec^, K =50 ft-lb and 0^=75 ft-lb-sec 
is plotted with only the indicated parameter varied, ceteris 
paribus . 

C. Bi\!.IPING AFREI) TO RECOVERY ROE 

In a further effort to reduce the maicimum swing, 0max, 
of the recovery system, damping wa,s s.lso ad.ded to the 
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VARIATION WITH Ci DAMPING 
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Figure 4 Effects of Damping Ci on 
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Figure 6 Effects of RPV Mass on with Cj Damping 
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Figure 7 Effects of Horizontal Member Length R- on 0j^3^ D^niping 
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Figure 8 Effects of Moment of Inertia on 0niax Cj Damping 



recovery rod. The added damping was included in the system's 
model, without cross coupling, by modifying the Rayleigh 
dissipative energy function to the form of: 

D=-g-C ^ 0 +-g-G 2 0 (35) 

The additional dissipation term precipitated through the 
equation of motion derivation described in section II, C. as 
the addition of a 0^0 term to the 0 eq.uation of motion, 
equation 23. This change was also easily incorporated into 

the computer program as indicated in the Appendix, 

An estimate of C 2 v/as obtained in a manner similar to 
that described in the previous section using the simple 
pendulum equation: 

0 +2 ^ + 00 ^0= 0 ( 26 ) 

where 




Only a slight amount of damping, 02=3.24 ft-lb-sec, 
added to the recovery rod of the base configuration, in 
addition to the horizontal member damping, was sufficient 

to reduce 0max to approximately 90 degrees at RPV approach 
speeds of 50 knots ( 0 q= 14 rad/sec). This can be seen by 
referring to Table III where 0 is listed as a function of 
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VARIATION WITH Ci, C2 DAMPING 
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•time after impact for approach, speeds ranging from 7 knots 
( 00=2 rad/sec) to 50 knots ( 0^=14 rad/sec). The maximum 



50-knot RPV approach speed v/as 50 degrees and 63 degrees, 
respectively. 

Further study of the modified base configuration v;as 
performed to obtain the magnitude of the RPV's linear 
acceleration. The absolute value of the RPV's linear 
acceleration was obtained using the following equation: 



where x, y and z are defined by equations 7a, b and c. A 
plot of equation 37 as a function of time for the modified 
base configuration is presented in Figure 9. It may be 
noted that the magnitude reduces to below one "g" within 
one second after the RPV contacts the recovery system. 

The sensitivity of 0max to each of the system para- 
meters is indica,ted in Figures 10 through 16, where the 
base configuration of R^=28 ft, R2=6 ft, il=4.66 slugs, 
Iq=1159 ft-lb=sec^, K=50 ft-lb, C^=75 ft-lb-sec and 0^= 

3.24 ft-lb-sec is plotted with only the indicated parameter 
varied, ceteris paribus. It can be noted from Figure 12 
that the system was sensitive to the recovery rod length, 

R 2 , v;hich may be interpreted as an implication that the 
engagement point for RPV capture on the vertical cables 
has considerable importance to system dynamic response. 



magnitude of 0 and for the modified configuration at the 




(37) 
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Figure 9 Linear Acceleration for Base Configuration 
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Figure 10 Effects of Damping C-, on 0j„ax with C2 Damping 
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Figure 11 Effects of Damping C 2 on 0niax 
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Figure 12 Effects of Rod Length R 2 on 0jjjgx C^, C 2 Damping 
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Figure 13 Effects of RPV Mass on with C, , C 2 Damping 
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Figure 14 Effects of Horizontal Member Length on with C^, C 2 Damping 
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Figure 15 Effects of Moment of Inertia Iq on 0 njax with Ci, C2 Damping 
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Figure 16 Effects of Spring Constant K on with C 2 Damping 



The system was also sensitive to the value of C^, but 
more sensitive to C2. This, however, should not present a 
problem, as and will be set by the design and will 
not vary in operation. 

The system was relatively insensitive to variations of 
the point along the horizontal member, R-j_, where the RPV 
was caught. The system was also relatively insensitive to 

f 

variations in the moment of inertia, I^, of the horizontal 
member and the mass, Li, of the RPV. The ma::imum displace- 
ment of 0 was not affected at all by changes in the spring 
constant, K. 

D. IivIPLEivISNTING- DAIIPING OR A RECOVERY CABLE 

The damping desired a,t the pivot of the horizontal 
member can be readily implemented in actual construction. 
However, it was not readily apparent hov/ the damping of 
the vertical rod, vvhich was a cable in the actual system, 
could be implemented. 

Two simple penduliim designs were analysed and modeled 
using the same techniques described to study the recovery 
system. The intent was to see if pivotal damping could 
be applied effectively to control the swing of a cable. 

The first design is shown in Figure 17 and. simulates 
a cable design for the vertical massless recovery rod that 

v/as used in the recovery system analysis. 7/ith no pivotal 
damping, an initial velocity v/as applied to to make 

*^2max ^ high angle w'ithout exceeding I80 degrees. 
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Then damping was applied to the pivot of the cable support 
to judge its effect on reducing the magnitude of 
The results are plotted in Figure 18. This approach of 
applying pivotal damping did not effectively reduce 
a,s desired. 

Another design, shown in Figure 19, v/as also analysed 
to see if it could, through pivotal damping, control the 
maximum magnitiide of 0 ^. Again the design simulated a 
cable instead of the massless recovery rod used in the 
analysis. A spring was used a.t the pivot to give a 

near horizontal static deflection. With no damping applied, 
an initial velocity was found for M 2 that caused 02 max 
approach I 80 degrees. Pivotal damping was then applied 
to 1,1^ to observe its effect on the magnitude of 02 * 
results shown in Figure 20 indicate that pivotal damping 
applied to this pendulum design was also ineffective in 
reducing 

V/hile neither of these designs demonstrated effective 
control of 0 ^, they by no means exiiaust the possible 
methods of implementing damping in a cable design, and 
further study is possible. 
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Figure 18 Results of Pendulum A Analysis 
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Figure 20 Results *of Pendulum B Analysis 



IV. CONCLUSIONS AND RECOL&iENDATIONS 



For the carousel recovery system to be feasible at 
approach speeds of 50 knots, a relatively small amount of 
damping v;as required at the pivot of the recovery rod in 
addition to damping at the pivot of the horizontal member. 
The damping at both pivots was required to prevent the 
maximum swing of the point mass RPV from exceeding approx- 
imately 90 degrees in the 0 direction, the coordinate that 
exhibited the largest change in magnitude. 

The recovery system v/as sensitive to the vertical 
height of the RPV at the time it was captured. The greater 
the distance, R 2 , between the ’RPV and the horizontal member, 
the greater the swing of the RPV relative to the horizontal 
member. V/hen selecting parameter values such as the damp- 
ing constants under the most severe conditions, the RPV 
should be assumed to have contacted the recovery system at 
the maximum length of the recovery rod or cables. 

The recovery system was relatively insensitive to 
variation in the distance the RPV was from the vertical 
mast when it was caught. The system was also relatively 
insensitive to variations in the mass of the RPV and the 
moment of inertia of the horizontal member. 

Variations in the spring constant for the torsional 
spring located at the pivot of the horizontal member had 
little effect on the swing of the RPV. Since it was not 
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important to controlling the motion of the RPV, it could be 
removed from the design. 

The magnitude of linear acceleration exhibited by the 
RPV, with the recovery system design that incorporated both 
horizontal member and recovery rod damping, was within 
Aquila load limits. 

Based on this analysis the carousel recovery system is 
still viable; however, it is apparent that the system will 
be sensitive to approach techniques used during the recovery. 

It is recommended that this analysis be continued, in 
order to treat the RPV as a distributed mass with its ovai 
moment of inertia and coordinate system, and caught by two 
recovery rods or. cables, one for each wing tip, instead of 
one. Expanding the scope of the recovery system model would 
provide a more detailed analysis of the RPV's linear lateral, 
longitudinal and normal accelerations and hence insure that 
they were individually within the RPV's design limits. 

It is also recommended that the effects of an RPV 

approaching at angles to the recovery system be investi- 
gated; this should result in increased motion in the ij/' 
direction. 

It is further recommended that continued study should 
be made of a means of applying damping to the swing of the 
vertica.l recovery cables. 
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APPEIIDIX 



I.IODELIKG PHE EQUATIONS OF I.IOTIOH USING- GSI.-IP 



An application-oriented soft\va,re package called the 
Continuous System Modeling Program (GSI.IP) v/as used to model 
and solve the nonlinear, cross-coupled equations of motion. 
GSI.'EP is a powerful softv;are package that can handle this 
type of problem routinely. In the analysis performed in this 
thesis two integrations using a fourth order variable inte- 
gration step Runge-Kutta technio^ue on each of the three 
equations of motion, Ref. 5. This sophisticated integration 
technique has the advantage of automatically adjusting the 
time increment of integration to meet the demands of the 
dynamic conditions of the simulation. The absolute vsXue of 
the estimated integration error and the relative m.agnitude 
of the estimated error are compared with user specified 

error-bounds. The step size is then adjusted to meet the 
desired error criteria. In the analysis performed in this 

thesis the default value of .0001 for both the absolute and 
relative errors was used. 

The program structure of GSI>IP is composed of three 
segments; Initial, Dynamic, and Terminal. Generally, data 

statements will appear in the Initial segment. In addition, 
calculations that are required to be performed only one 
time during the simulation can be conveniently placed in 
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this segment. The Dynamic segment is usually composed of 
structural statements that describe the dynamics of the 
system or explicitly describe a set of differential equa- 
tions. The Terminal segment is the last segment in the 
program and is usually made up of output control statements. 
The actual program used in this analysis for the un- 
damped system will now be described. List of symbols used 
in the Computer Programs; 

R1 - Length of horizontal member in feet 

R2 - Length of recovery rod in feet 

M - Mass of RPV in slugs 

HIO - Moment of inertia of horizontal member in ft-lb- 
2 

sec about its pivot point 

K - Spring constant, in f t-lb/rad , • of spring at pivot 

of horizontal member 

X5I - Initial angular velocit 3 /‘ of the recovery rod in 
the 0 direction in units of rad/sec also equal 
to approach velocity of RPV in ft/sec divided 
by R2 

G- - Acceleration of Gravity 32.2 ft/sec^ 

XI - 0 

X2 - 0 

X3 - 1/^ 

X4 - 0 

X5 - 0 

X6 - ijr 

X4D0T - 0 
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X5D0T 


- 0 


X6D0T 


- f 


CO 


- G^, damping constant in ft-lb-sec for the damper 
located at the pivot of the horizontal member 


C02 


- C^, damping constant in ft-lb-sec for the damper 
located where the recovery rod attaches to the 
horizontal member 


AD 


- Determinant of A as defined in Section II. C. 4 


ADI 


- Determinant of A1 as defined in Section II. C. 4 


AD2 


- Determinant of A2 as defined in Section II. C. 4 


AD3 


- Determinant of A3 as defined in Section II. C. 4 


B1 


- is b.^ defined in Section II. G. 4 


B2 


- is b 2 defined in Section II. G. 4 


' B3 


- is b^ defined in Section II. C. 4 


A1 


- is a^^ defined in Section II. C. 4 


A2 


- is a ^ 2 defined in Section II. C. 4 


A3 


- is a ^2 defined in Section II. C. 4 


A4 


- is defined in Section II. C. 4 



A, B, C, D, E, P, H, L, S, T - v/ere used to simplify 

program notation and are 
defined in the program 

A description of the program shovni in Figure 21, follows: 
Lines 1-3: llalce up the initial segment, para-meters 

that 'would remain constant during the 
integra,tion v;ere asGigned here. 

Lines 4-34: Hake up the Dynamic segment. 

Lines 5-13: Simplifying notation. 
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Figure 21 Computer Program 



Lines 14-28: 



Lines 29-31J 

Lines 32-34: 

Lines 35-41: 
Line 36: 

Lines 38-39: 



The algebra used to solve for 0 , 0 and Ip' 
prior to integration as described in 
Section II, C. 4. 

Perforin the first integration from angmilar 
accelerations to angular velocities. 

Perform the second integration from angular 
velocities to angular position. 

Make up the Terminal segment. 

Limits the time period that will be modeled. 

Prints output as a function of time. 



When damping v/as added to the horizontal member the 
only changes that occurred to the program v/ere the follo’vving 



Line 2: 


add 


LT\ 

!! 

O 

O 


Line 23: 


add 


-CO * X4- 


When d. 


amping 


was also added to the recovery rod the 


additional 


changes that occurred to the program were as 


follows : 






Line 2: 


add 


,C02 = 3.24 


Line 24: 


add 


-C02 * X5 
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